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Development of quantum engineering put forward new theoretical problems. Behavior 
of a single mesoscopic cell (device) we may usually describe by equations of quantum 
mechanics. However if experimentators gather hundreds of thousands of similar cells 
there arises some artificial medium that one already needs to describe by means of new 
electromagnetic equations. The same problem arises when we try to describe e.g. a 
sublattice structure of such complex substances like perovskites. It is demonstrated that 
the inherent primacy of vector potential in quantum systems leads to a generalization 
of the equations of electromagnetism by introducing in them toroid polarizations. To 
derive the equations of motion the Lagrangian and the Hamiltonian formalisms are used. 
Some examples where electromagnetic properties of molecules are described by the toroid 
moment are pointed out. 



1 Introduction 



Let us remind a known thing that says "it is impossible to introduce electrodynamics of 
matter in general" (E. A. Turov, 1983). For example, different types of crystalline 
structures of matter lead to the alignment of one or other type of polarizations in the 
matter considered. So the necessity to introduce in the equations of high tensor polar- 
ization is arose. But the most intricate case is the toroid polarization one. The fact is 
that the ideal static toroid moments do not interact with each other. So the dynamic 
alignment of toroid moments is impossible thanks to electric and magnetic interactions. 
But this takes place for example in perovskites [1, 2] and have to be explained. The other 
case when the local toroid moments can align is connected with their proper oscillations 
that permit them interact with each other [3] through electric and magnetic fields. 

It is noteworthy to remark that toroid moments are the multipole sources of field-free 
vector potentials. Therefore, electromagnetotoroidic equations we are forced to express 
in terms of vector potentials. The first part of our report is devoted to introduction of 
the equations mentioned in Lagrangian formalism. In the second part we introduce toroid 
moments in Schrodinger equations. 

2 Toroid moments in Lagrangian formalism 

2.1 Canonical formulation of electrodynamics 

We begin with usual canonical formulation of electrodynamics (e.g. [4]). The inter- 
acting system of electromagnetic field and non-relativistic charged particles is specified 
partly by a discrete set of variables, namely the coordinates of the charged particles, 
and partly by a continuous set, which we take to be the values of the vector potential 
in the Coulomb gauge. The Lagrangian L will thus be a functional of A and A if the 
particle coordinates and their velocities are fixed, and a function of the q a and q a if the 
vector potential and its time derivative are fixed. We write L = L[A, A; q, q}. In the 
application of Hamilton's principle, the particle and the field coordinates are to be varied 
independently. The Lagrangian must then be chosen so that variation with respect to the 
particle coordinates gives Newton's law 



and variation with respect to the field coordinates (subject to the Coulomb gauge condi- 
tion) gives the equation of motion 



m a q, 



ioc = e Q (E + -q a x B) 



V 2 A - Ia = - 
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for the vector potential. 

A suitable Lagrangian is obtained by setting (see e.g. [4]) 
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- (curlA) 2 ]^ (2-3) 
- / J(r)-A(r)dV = J2-^-MQa,t). (2.4) 



Here L par is the Lagrangian appropriate to a system of charged particles interacting solely 
through instantaneous Coulomb force; it has the simple form of "kinetic energy minus 
potential energy" . L ra( i is the Lagrangian for a radiation field far removed from the charges 
and currents, and has the form of "electric field energy minus magnetic field energy". The 
interaction Lagrangian L int couples the particle variables to the field variables. The part 
of the Lagrangian that involves the vector potential can be written as the integral over a 
Lagrangian density C, 

1 A 2 1 
£=-[—- (curl A) 2 ] + -J^rj-AW (2.5) 

07T & C 

Since the Coulomb gauge condition is being used as a constraint, the transverse current 
density J -1 has been substituted for the total current density J without affecting the 
Lagrangian. 

To verify that L given by Equation (2.1) is indeed a suitable Lagrangian, we write 
down the Euler-Lagrange equations, beginning with those for the particle coordinates. 

d PL _ dL_ 

dt dq ai dq ai 

In this particular case we have 

m a q ai - e a [j2 , H |3 (g« - QfSi) ~ \ aj - ~ lr~) = °" (2 ' 6) 

^ a l^a Q/jI c j c x oq ai oq a j' 

Using the expressions B = curl A and E" = — grad0, E^ = — ^A with 0(r, t) = 
J2 a \q Ca -r\ f° r tne transverse fields in terms of the vector potentials, we find that 

m a q a = e a E{q a , t) + ^q a x B(q a , t) (2.7) 

which is nothing but the second law of Newton with the Lorentz force. 

To obtain the field equations, the functional derivatives of L must first be calculated 
from the Lagrangian density C. The Euler-Lagrange field equations 

d 5L 5L 

in this case give 

curl B = - E + — J 1 - 

c c 
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which leads to the usual evolution equation for the vector potential in the Coulomb gauge: 



curlcurlA + \ A = —J 1 . (2.8) 

cr c 

The conjugate momenta corresponding to the Lagrangian (2.1) are defined in the usual 

way as 

Pa = w = m< ^ + v A(q ' t} (2 - 9) 

and 

nW = t^TT = ifi(= -E(r)/4vrc). (2.10) 



9A(r) 4ttc : 
Proceeding in the conventional way we get 



H[U,A;p,q] = £ p a • q a + / II • ArfV - L 



E^-[P Q --A(q,t)] 2 + iE ? 



a 2m a Q c 2 ^ |q a - q^| 

+ — /[(47rcn) 2 + (curlA) 2 ]rf^ (2-11) 

The Hamiltonian equations reproduce the equations of motions (2.7) and (2.8) as they 
must be. The pair of equations 

dH ■ SH_ 

dp ai dUi 

are equivalent to the relations between the canonical momenta and Lagrangian velocities, 
where as the equations 

dH ■ 5H 

dq ai bAi 

yield once again the Lorentz force law (2.7) and the inhomogeneous wave equation (2.8) 
for the vector potential with the transverse current density as source. 



2.2 Formation of an equivalent Lagrangian 

It is well known that in classical dynamics the addition of a total time derivative to 
a Lagrangian leads to a new Lagrangian with the equations of motion unaltered. La- 
grangians obtained in this manner are said to be equivalent. In general, the Hamiltonians 
following from the equivalent Lagrangians are different. Even the relationship between 
the conjugate and the kinetic momenta may be changed. 

More over, let us notice that the basic equations of any new theory can not be intro- 
duced strictly deductively. Usually, either they are postulated in differential form based 
on the partial integral conservation laws or transformations of basic dynamical variables, 
whose initial definitions usually have some analogs in mechanics. We will mark that we 
need to do so not only by inertia of thinking but also because of the fact that most of 
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our measurements have its objects as individual particles or use them as test one. The 
situation is the same in electromagnetism and in gravitation. Attaching of geometri- 
cal meaning to any dynamical variables plays the crucial role in this case. Even that 
slight modifications of electrodynamic equations that we performed, demands the clear 
understanding of this situation. 

As a simple most example we will consider the "geometrization" of momenta in non- 
relativistic quantum theory when the external classical electromagnetic field, acting on 
free particle, is available. Let us begin with the classical Lagrangian: 

£(r,v;t) = l -mw 2 - g[0(r;f) - v • A(r;*)], 

that leads to the equations of motion of particle with charge q and mass to, moving under 
the influence of electric and magnetic potentials and A. To move up to Hamiltonian 
formalism, more closer to the quantum one, let us find the canonical momenta of the 
particle: 

p := 8 V jC(t, v; t) = tov + qA(r; t). (*) 
The Hamiltonain in this case is defined as: 

H := pv - £ = ^-[p - qA] 2 + q< j>. (**) 

Let us use the initially defined mechanical momenta 77 := tov (gauge invariant quantity 
that can be directly measured!) to write the expression in square bracket in (**). We find 

[p - gA] 2 = [mv + gA-gA] 2 = 77 2 . 

It is of course the simplest tautology, nevertheless under introduction of electromagnetic 
field in Schrodinger equation it leads to some subtlety. Quantum physicists random say or 
write that the gauge invariant introduction of external electromagnetic field in Schrodinger 
equation is occurred by extended momenta (of particle), implying the substitution p — > 
P - ?A. 

In fact the operation to move up from the classical description to the quantum one 
with the introduction of electromagnetic potential we imply the following substitutions: 

77 := tov = V to p = 77 + qA — > V 

and 

. In this the mathematical expectation of meansquare deviation between "geomtrized" 
momenta (operator of spatial displacement of particle) and field (potential) angular mo- 
ment of the particle gA are equated to the quantum mechanical expectation of evolution 
operator of the particle considered: 

d 1 

* = — * + [(mV + gA) 2 + qd)W. (***) 
at 2m 
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In deduction of Schrodinger equation in [5] the similar view is taken into account. Power 
a.o. [5] (see also [14]) begin the introduction of electrical polarization in Schrodinger 
equation starting with the equivalent Lagrangian. 
An equivalent Lagrangian to that of (2.1) is [5] 

L new = L - ~[ / p ( r ) • A ( r ) dV \- ( 2 - 12 ) 
cat J 

Taking into account the Coulomb gauge condition divA = and toroid contribution we 
may substitute 

P =>• P x + P 11 + curlT 6 , (2.13) 

with the contribution of P" being vanished. 

Thus we have the following equivalent Lagrangian 

L new = L - -^-[f (P ± (r) + curlT e (r)) • A(r) dV]. (2.14) 
Writing it in the explicit form we get 

L new = L par + L rad + - c J (J x (r) - P X (r) - curlt e (r)) • A(r) dV - 

- - J (P x (r) + curlT e (r)) • A(r) dV (2.15) 

The field conjugate to the vector potential A is now 

n(T) _ dLneW _ E± W + 4vr(P ± (r) + curlT e (r)) 

1 ' dk(v) 4vrc 1 ' 1 

In fact, if we define the new vector field in a medium D by the relation 

X>(r) = E ± (r) + 4 7 r(P ± (r) + curlT e (r)) = D x (r) + 47rcurlT e (r) (2.17) 
for the conjugate field we get 

n(r) = -^ (2.18) 
Aire 

Now from the definition of T> we get 

curlX>(r) = --B(r) + 4tt (curlP(r) + curlcurlT e (r)), (2.19) 

c 

under curlE(r) = — - c B(r), since only the free field E is generated due to the change of 
magnetic field B as before. 

The new Lagrangian is a function of the variables q Q , q a and a functional of the field 
variables A, A, and the equations of motion follow from the variational principle. We 
have 

ajnew 1 i 

(J x (r) - P (r) - curlt e (r)) = - (ccurlM ± (r) - curlt^r)), (2.20) 



<9A(r) c ' c 



where the relation ccurlM ± (r) = J" L (r) — P (r) is taken into account and 

Qjjnew 

WMrjJd^) = A^ £l]kBk (2 ' 21) 
as usual. Applying the Euler-Lagrange equations of motion 

d dL new ^ d dL new dL new 

dt dljr) + yd^j didAitf/dxj) ~ 3A~(r) = ° (2 ' 22) 

we obtain 

- ?iEl + -LcurlB(r) - - (ccurlM ± (r) - curlt e (r)) = 0. (2.23) 
4tt c An c 

Introducing the following re-notations, corresponding to the initial identification by Maxwell 
for local conastrain between the derivatives of electric and magnetic fields: 



c 

we may rewrite (2.23) as 

1 



-T e ' m !„=}► T curlT m ' e \ n (2.24) 



curlB(r) = -2>(r) + 47r(curlM ± (r) + curl curlT m (r)) (2.25) 

The relation (2.24) demands some comments. Both T 6 and T m represent themselves in 
essence closed isolated lines of electric and magnetic fields. So they have to obey the usual 
differential relations similar to the free Maxwell equations (see in [2] and [6]). However, 
remark that signs in (2.24) are inverse in comparison with the corresponding Maxwell 
equations because the direction of electric dipole is accepted to be chosen opposite to its 
inner electric field [7]. 

If we define the auxiliary field H to be 

H(r) = B(r) - 47r(M ± (r) + curlT m (r)), (2.26) 

then it deduces to 

curlH = -t>. 

c 

But the latter formula is unsatisfactory from the physical point of view. It is easy to 
image the situation when B and M -1 are absent, because the medium may be composed 
from isolated aligned dipoles T m [1, 2] and each T m is the source of free-field (transverse- 
longitudinal) potential but not B [3]. So the transition to the description by means of 
potentials is inevitable. 
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2.3 Hamiltonian for new Lagrangian 

Proceeding as before we get 

H new [U,A;p,q] = ]Tp Q -q a + [ U ■ AdV - V 



t 2m Q c ' 2 ^ |q a - 

+ 77- / {[47r(P + curlT e ) - £>] 2 + (curlA) 2 ]^ 

87T J 

+ i y (P X + curlt 6 ) - AdV". (2.27) 

Taking into account ^P" 1 = Ij^ - curlM" 1 and the relation (2.24) and also div[M x A] = 
A • curl M + M • curl A the fourth term of the foregoing equality can be reconstructed 
after integration and H can be presented as 

H^[n,A;p, g ] = E^-[P a --A(q,i)f + ^E ' 



« 2m * a c " v ' 2 loa - q^l 

+ /{[47r(P + curlT e ) - T>} 2 + (curl A) 2 } rfV 

57T J 

+ i y J x - AdV - JlSA-BdV - yB-curlT"W. (2.28) 
The Hamiltonian equations in this case give 

?« = |^ = E — (P« - — A(q Q , 0), A = ^ = c[4tt(P + curl T e ) - T>\ 
which together with 

n, = - ^ = E — (P« " — A(Qa,*)) - -(P + curlt 6 ), - J-[V x (V x A)],, 
oAj ^ cm a j c c 47r 

dfe aCm ai Pm c 1 qa ' <9g a / 2V|qa-q/3| ga QP 

once again yield the Lorentz force law (2.7) and the inhomogeneous wave equation (2.8) 
for the vector potential. 



2.4 Canonical quantization 

Let us now draw the Schrodinger picture operators corresponding to the quantities ob- 
tained above. To do so we need to interpret the coordinates q a of the particles together 
with their conjugate momenta p Q and the coordinates A(r) of the field together with 
their conjugate momenta II(r) as operators in Hilbert space. Thus we have six opera- 
tors q ai and p ai , (i = 1,2,3) associated with each particle a and six operators A^r) and 
Ilj(r), (i = 1,2,3) associated with each field point r. These operators are Schrodinger 
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picture operators since they are time-independent. The corresponding Heisenberg picture 
operators may be written as q a (t), p a (t), A(r, t), II(r, t) respectively. 

The Schrodinger operators are to satisfy the following canonical commutation rela- 
tions, precisely 

[qai,qpj] = 0, \pd,ppj] = 

[QahPfSj] = ih5 afS 5ij (2.29) 

for the particle variables and 

[Mr), Mr')} = o, [n^n.y)] = 

[^(r),H,(r')] = ift^(r-r') (2.30) 

for the field variables. These equations are to hold for all a and (3 and all % and j. The 
commutator bracket of any particle variable with any field variable is to vanish. The 
canonical commutation relations are also satisfied by the Heisenberg dynamical variables, 
provided all operators are evaluated at the same time t. These follows from the equations 
(2.29) and (2.30) and the relation between two pictures: 

n(t) = e m)Ht ne -m)m ^_ 31 ^ 

where H is the hamiltonian operator and Q(t) is the Heisenberg operator corresponding 
to the Schrodinger operator Q. The operator Q(t) satisfies the Heisenberg equation 

Q(t) = ^-[Q{t),H], (2.32) 
in 

while the Schrodinger picture operator corresponding to Q(t) is given by 

Cl = 4 \n,H}. (2.33) 
in 

Taking the expression (2.28) for the Hamiltonian one can easily verify that the equation 
(2.32) together with the commutation relations leads to the expected equations of motion. 



2.5 Formal deduction of equations of electromagnetotoroidics. 

Let us suppose that in an electromagnetic medium there are no free charges and 
currents. So we may rewrite usual Maxwell equation in the following symmetrical form: 

1 • 

curlB = -D + 47rcurlM, B = H + 4irM in the whole 5r, 

c 

curlD = -^B + 4vrcurlP, D = E + 4ttP in the whole 3? 3 , (2.34) 

being of these equations interchange to each other as before through the self-reciprocal 
exchanges B — > =fD, D — > ±B, M — > =pP, P — > ±M and the conditions divB = 
and div D = are fulfilled. Hence is valid the following changes [8] 

B =^ curl a m + a e , D =^ curl a e - a m , 
curl M =>> curl M + curl curl T m curl P =>• curl P + curl curl T e (2.35) 
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whereas the last substitution has already been implied canonically in the fore going para- 
graphs. As a result we obtain for basic equations 



curl curia" 1 + a m = 47r(curlM + curl curl T m ), 

curl curia 6 + a 6 = 47r(curlP + curlcurlT 6 ), (2.36) 

It is necessary to emphasize that the potential descriptions electrotoroidic and magne- 
totoroidic media are completely separated. The properties of the magnetic and electric 
potentials a m and a 6 under the temporal and spatial inversions are opposite [2] . We also 
see that the substitutions (2.28) in (2.25) produce (2.34) as well. 

If divD ^ and in the medium there does exist free current we have to remember 
Dirac's approach to the constrained systems and use along for instance [9] its application 
to electrodynamics of continuous media. 

2.6 Solution to the electromagnetotoroidic equations 

Let us first solve the static equations. In this case we write 

curlcurl(a m - 47rT m ) =: curlcurla" 1 = 4 7 rcurlM ± (r), (2.37) 

curlcurl(a 6 - 4ttT 6 ) =: curl curia 6 = 4vrcurlP ± (r), (2.38) 
here a is an analog of H. From the vector analysis we can see 

curlcurla" 1 = grad diva m - V 2 a m = -4tt grad divT m - V 2 a m (2.39) 

where we use the condition diva™ = The solutions of the equations (2.37) and (2.38) 
can be written as 



a m 



r curl'M(r / ) , . r div'T m (r'j , . 

/ -i -7T- d V + V / dV 2.40 

J r — r J r — r 



Analogically we get 



a e 



J r — r J r — r 



Now let us return to the general case. Adding — (Aix/c 2 ) t" 1 ' 6 in the both sides of the 
equations (2.35) and (2.36) for a we get the wave equations 

curl curia™ + 1 a m = 4 7 rcurlM ± (r) - -^T m (r) (2.42) 

c c 

and 

curl curia 6 + \ a 6 = 4vrcurlP ± (r) - -^T e (r) (2.43) 

c c 

The solutions to these equations can be written as [10] 

[-V'(V • T m/ ) - 4tt V x M') + (47r/c 2 )T' 

all space 



a m 



= _J_ f hVjV • Tm/ ) ~ 4vr V- x M') + (47r/c 2 )f m/ ] dv , 
Att J |r - r'l 
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in. a 



47T 7 |r - r'| 

aii space 

here M' denotes M(r'). Under the Coulomb gauge the substitution div a m,e = — 47rdivT 
is valid. The brackets [] in the above equations are the "retardation symbol". This sym- 
bol indicates a special space and time dependence of the quantities to which it is applied 
and is defined by the identity 

[/] = f(x',y',z',t-r/c). 

It is obvious from the foregoing solutions that the quantities that an experimentator 
measures are at and a and the fields, they generate, differ from the old ones. It leads to 
the fact that we no longer work with B and D but with some new quantities which need 
to be defined differently. So instead of (2.35) we should now write 

(3 = curia™ + a e , 6 = curia 6 - a m . (2.46) 

Note that the quantities (3 and S may be measured and possess physical significance. 
For example, the experimentator, measuring magnetic field in a needle-like hole, dug 
in magnetic, not along the principal axis of magnetization, beside the contribution of 
Ampere current of magnetization curlM measures the contribution of inhomogeneous 
toroidazation that is obvious from (2.44) and (2.45). Thus introduction of potential and 
magnetic field based on the magnetization of media only, generally speaking is incorrect 
[3]. 



3 Multipolar interactions in Schrodinger equations 



A large number of works have been devoted to this problem. Unfortunately most of 
them are confusing. We demonstrate it on the basis of detail paper by K. Haller and 
R. B. Sohn [11]. The matter is that they begin with the expression independent of the 
scalar- longitudinal contributions (p, divj) and after the integrations by parts in the final 
expression restore them. It is a very common methodical error. 

3.1 Description of non-relativist ic interaction of electrons and 
photons 

Bere we proceed and set the problem in terms used in [11]. The Bamiltonian that describes 
the interaction between photons and non-relativistic Schrodinger electrons is given by 

H c = H - J J(r) • A T (r) dr + ±- f p{v) A T (r) • A T (r) dr + 

+ / drdr', (3.1) 

J 8tt r - r' ' K J 
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where A T is the transverse vector potential with div A T = 0. H is the Hamiltonian 
for noninteracting electrons and photons that can be represented by 

H = H (e) + #0(7), (3-2) 

where 

H (e) = J ^+(r)[-(2m)- 1 V 2 + V(r)]^(r)rfr (3.3) 

with V(r) representing an external short-range potential (for example, the shielded Coulomb 
potential of a static nucleus). H (j) is the Hamiltonian for free transverse photons and 
is given by 

#o(7) = \j [E T (rf + B(r) 2 ]dr, (3.4) 

where E T (r) and B(r) represent the transverse electric and magnetic field, respectively. In 
equation (3.1) p(r) represents the charge density p(r) = e ?/> + (r)-?/>(r) and J(r) represents 
a current 

J(r) = — [V^ + (r)^(r) - ^+(r) V^(r)]. (3.5) 
The commutation rules of these fields are 

{^ + (r),^(r')} = 5(t - r') (3.6) 

for the electron field; for the photon field, — A T (r) has the transverse electric field E T (r) 
as its conjugate momentum, so that the nonlocal commutation rule for the transverse 
components is given by 

WWW - -CMC - O + Irl-i^)- 

The current J(r) is conserved under the time dependence provided by the Hamiltonian 
so that 

divJ(r) = -i[H ,p(r)]. (3.8) 

The current 

j(r) = J(r) - (e/m)p(r)A T (r) (3.9) 
is conserved under the time dependence provided by the Hamiltonian so that 

divj(r) = -i[H c ,p{v)\. (3.10) 

H c may be expressed as 

H c = H - J [J(r) + j(r)] • A T (r)dr + / dvdv'. (3.11) 

When ip is quantized in orbitals that are solutions of 

[-(2m)- 1 V 2 + V(r) - u^Ui = 0, (3.12) 
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(where C/j includes both bound and continuum states), H Q (e) is given by 

H o( e ) = e n e nU n , (3.13) 

n 

where e+ and e n designate creation and annihilation operators, respectively, for electrons 
in the n orbitals. #0(7) can be represented as 

#0(7) = E / dk|k|aj i) (k)a e(i) (k), (3.14) 

i=l,2 J 

and describes free photons in the two polarization modes. 
3.2 Introduction of toroid moments 

A formal procedure was first proposed by E. G. P. Rowe [12] permits to obtain the 
complete solution to the problem of multipole expansion of electromagnetic current [3], 
which consists of the replacing of some vector function of current by the three (in general 
unlimited) series of multipole parameters. The multipole parametrization of interaction 
Hamiltonian of an arbitrary system with external fields under the Coulomb gauge has the 
form [3] 



(21 + 1)!! ./2Z + 1 



00 1 00 

Hc = ~SJl,£2 B »K2i + 2n + l)!!V 4vr 

x {r 1 M^ n \t)Y lm (W)A n (r-B)\ r=0 + (3.15) 
+ U[Qim{t)8 nfl l\- 1 - T%(t)]Y lm (V)A n [(l/c)r • D + (47r/c)r • j] r=0 }. 

where Qi m (t), connected with Coulomb multipole moments of the charge distribution of 
the system are 

Q lm (t) = V^d J r^tU^jM)^, (3.16) 
(t) are the magnetic multipole moments or their radii 



«S°W = T V 7TT V 2(TT / ^mr,t)1r, (3.17) 

and T^\t) are third family of multipole moments [13], the toroid moments and their 
radii, namely 



T^\t) = -^fyj / r^frUJr) + 2 -^^Y* +lrn (r)] ■ j(r,f) dr, (3.18) 

In the case considered by Haller and Sohn [11] we have to substitute into (3.15) 

Q lm (t) = ^ = i[H c ,Q lm ] (3.19) 
at 
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and 

j = J - epA (3.20) 

in the formulas (3.15), (3.16), (3.17) and (3.18). In this approach there does not appear 
any longitudinal contributions that are fictitious (mutually cancelled out) in the expres- 
sion, deduced by Haller and Sohn. The formulas (3.16 - 3.18) obtained by us give correct 
multipole parametrization for the interaction energy transverse degrees of freedom of a 
non-relativistic system with radiation field. 

3.3 Toroid moments in Schrodinger equation 




Figure 1 



Now suppose that we 
consider some electron 
involved in the molecu- 
lar structure. We may in- 
troduce a coordinate sys- 
tem as in Fig. 1. Then 
it may be shown that 
the Schrodinger equation 
for this electron interact- 
ing with external sources 
of electromagnetic fields 
has the following form 
[14] 



ity(q) = [-^ (q) f + V(q) + e 2 j 



(q')0(qO rfq 



q - q' 



- J P(r, q) • E ± (r) dr - J M(r, q) • B ± (r) dr + 



+ 



2^ [ I n(r ' q) 



xB(r)rfr] 2 0(q), 



(3.21) 



with n(r, q) = -f(q - R) 5(r - R). 

To get the Schrodinger equation avobe authors used the following substitutions in 
the Lagrangian for an electron field uncoupled from the radiation field: V is replaced 
by V + jp-A(q) and the wave function describing electron field initially is transformed 
according to 

^(q) = e- S( ^(q), 

where S(q) is given by 

5(q) = ±J P(r,q).A ± (v)dv. 
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There are two ways to introduce here the toroid contributions. The first one is straight- 
forward to use the substitution (2.30) Then we obtain 

J P-E^dr =>■ J P-E ± rfr + J curlT e -E ± dr = J P-E ± rfr + J T e -curlE- L dr. (3.22) 

Analogically 

| M-Bdr^ | M-Bdr+ J T m ■ curlB dr. (3.23) 

More reasonable approach is developed in [2]. It goes back to the classic multipolar 
description of quasimolecular structure [9,15]. According to it we may use immediately 
the multipole expansion of the densities P(r, q) and M(r, q) as follows [2]: 



= - J d(r, q) • E ± (r) dr = -Q • E x - T e curlE 



-P e • curlcurlE - ^Q^V^f + V,E^) - 7^.V i (curlE- L ) i (3.24) 

where Q = / d(r) dr is the conventional total electric dipole moment of the system, 
T e = \ J r x d(r) dr is the axial toroid moment and (see [14] and Appendix) 



P(r, q) = -e{(q - R) - l(q x R) - 
[(q - R) t (q - R)j - \{q - R) 2 5 l3 ] + ■ ■ ■} S(r - R). (3.25) 



Analogically for M(r, q). 

Let us now define the axial toroid moment. From the classical point of view on the origin 
of electric dipoles we see that the multipole expansion of the distribution density P(r, t) 
contents generally three kinds of dipole moments (see Formula (3.24)). Working within 
the scope of classical framework after Power a. o. [5] we see that T 6 = for a separate 
"atom". Really, all the electric dipoles are characterized by its space vectors, with the 
origins lying in the origin of the coordinate system (center mass of the system considered) 
and the endpoints in the electron coordinates. The sum of these vectors characterizes the 
total electric dipole moment of the system. To demonostrate the forementioned conclusion 
we consider the following example. 
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Let us have a system con- 
taining N atoms where 
each j-atom itself forms a 
subsystem [see Fig. 2]. 
Let "O" be the center- 
mass of the system as 
a whole. If the system 
contains only one atom, 
say "A", which on his 
part consists of one nu- 
clei and n electrons, the 
electric dipole moment of 
this system relating to 
the point "O" can be 
written as 

Q = -e£(R- qi ), d = Q5(r-R). (3.26) 

i=l 

Then the axial toroid moment of the such a system with respect to its proper center-mass 
is 

T M = -J dxrdr = - £(R - qj x (R - Qi ) = 0, (3.27) 

as Ti = R — qj and R is fixed. If the molecule contains N atoms then we define the total 
toroid moment of the molecule as a whole [16] with respect to its center-mass by the 
following formula 

N n N n 

T (d) = ? E E(R* " %') x R . = 9 E Et% x R,"] ^ 0, (3.28) 
z j=l 1=1 z j=l 1=1 

although the axial toroid moment for each separate j-atom with respect to its proper 
center-mass is equal to zero as before. 

As a particular example where electromagnetic properties of molecules are described by 
the axial toroid moment, we point out the phenomenon of " aromagnetism" [17]. There 
exists substances with a closed chains of atoms like benzol ring. In the experiment, 
the microcrystals of the aromatic series substances (antrazen, fenantren etc.), suspended 
either in water or in other liquids, were reoriented by an alternating magnetic field so that 
the modulation of polarized light, propagating through the media given, was observed. 
This reorientations can be easily explained if the aromatic molecules are considered to 
possess axial toroid moment (see Fig. 3) as an electromagnetic order parameter. 
The latter can be clarified as follows. All molecular wave functions \P of the benzol ring, 
being the main fragment of aromatic molecules, are first classified through the irreducible 
representations of the point group symmetry and then the asymmetric representations 
(E 2g and E lu ) are selected among them. When such molecules are packed into a molecular 
crystal, the parallel orientation of their axial toroid moments are preferred and the crystal 
as a whole can possess a macroscopic axial toroid moment [18]. 




Figure 2 
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Figure 3. Here is a fragment of an aromatic molecule where Ci is a carbon atom, 
is a hydrogen atom or its substitute with C; L . The vectors dj denotes the electric dipole 
moments, being due to the molecular orbital with the symmetry E 2g of the representation 
of the point group D^. The molecular wave functions ^E 2g and ^E lu with this symmetry 
are built [18]. The molecules in this states are shown to have the axial toroid moment. 

Now let us see what " aromagnetism" is? Microscopic crystals of aromatics like anthracene, 
penthacene, phenanthrene and others were suspended in water or some other liquids. 
When alternating or rotating magnetic field was applied to the suspension the change 
of orientation of microcrystals was observed. This new magnetic property of aromatic 
substances was named "aromagnetism" [17]. This phenomena can not be explained by 
the standard way like ferromagnetism since the organic molecules do not possess magnetic 
moments neither orbital nor spin origin. 

The theory of aromagnetism origin was proposed in [18]. It was assumed there that 
the sample of aromagnetic substance is possessed of an axial toroid moment G (ATM). 
This moment describes the vortex distribution of the electrical dipoles d a deposited at 
points r a and equal to 

G = (l/2)5> a dJ, 

where the sum is taken over the whole sample. Electrodynamics theory of ATM was 
developed in [19]. The energy of ATM in external field is equal to 

<9TT 

£ = -GcurlE = G— . 

dct 

Thus toroid moment G, electrical by its nature, can interact with alternating magnetic 
field H(£). In the following we consider as 'aromagnetic' any molecule or macroscopic 
sample that possess toroid moment G. The main feature that singles out aro magnetics 
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among the other magnetics is their interaction with magnetic field time derivative dH/dt 
(not with magnetic field H itself). The theory previously given in [18] may be extended 
so that new substances with aromagnetic properties can be predicted. 

Indeed it was shown (see [18]) that the cause of molecular aromagnetism is the ex- 
istence of quantum molecular electronic states (MO) where the quantum mean value of 
toroid moment < G > of molecule is not equal to zero. In [18] only molecules with group 
symmetry were considered. It is obviously that molecules with the other symmetry 
groups may be examined. The scheme of the examination is the following. Toroid moment 
is a pseudovector and for any given group we know the irreps of its transformation under 
the group operations. For example in the group D ih the irreps of pseudovector are : A 2g 
for G z and E g for G x and G y components of G. Considering direct products of this irreps 
(A 2g and E g in a given case) with the other irreps of this group it is possible to find out 
the symmetry of MO where the mean value of toroid moment is not zero (E g or E u for 
group D±h). Then it is possible to find these MO in explicit form and calculate the mean 
value < G >, as it was described in [18]. Toroid polarization coefficients of the molecule 
can also be calculated in analogous way. 

Figure 4- Distribution of dipole moments of 
oxygen atoms of the molecule phloroglycine. 
It can be estimated from this figure that 
molecule has axial toroid moment. 

The other important case of molecular aromagnetism is the vortex distribution of 
dipole moments of oxygen or nitrogen atoms they have due to their lone pairs of electrons. 
On the Figure 4 the stereomer of the molecule of phloroglycine is represented and dipole 
moments of the oxygen atoms are shown. The stability of this stereomer was demonstrated 
by numerical calculation by standard molecular dynamic method. It follows from the 
picture that vortex distribution of the electrical dipole moments on the oxygen atoms 
exists and that this stereomer of phloroglycine is aromagnetic one. The other examples 
of organic molecules with the same kind of aromagnetism are given here. 

To explain structural aromagnetism let us consider molecular crystal of indandion-1,3 
pyridine-betaine (IPB) polar intramolecular salt (Figure 5, left). Due to charge transfer 
this molecule has electrical dipole moment. The fundamental unit of the crystal of IPB 
is shown in Figure 5 (to draw the Figure 5 the picture 1.9 from the book [20] is used) 
and on the figure we show the dipole moments of the molecules. It can be seen from the 
figure that toroid moment of the crystal cell is not equal to zero, i.e. this substance has 
aromagnetic properties which we call structural aromagnetism. 
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Figure 5. Structure of IPB molecule (left). Deposition of IPB molecules in fundamen- 
tal crystal unit (right). The electric dipole moments of the IPB molecules are shown that 
demonstrates it vortex (toroid) distribution. 

Let us see one more example. In a cubic perovskite the crystallographic plane (1, 1, 1) 
contains triangular sublattice of oxygen atoms [21,- 23] (see Figure 6). Magnetic moments 
of these atoms can be in two states of orientating symmetry that are separated by tem- 
peratural phase transition. 




Figure 6. Structure of type perovskite-type 
compound Barium Titanate BaTiO^- In re- 
ality atoms fill up the major part of space. 
Here the position of their center is shown 
only. 
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Figure la. Magnetic moments of oxygen atoms are orientated in such a way that they 
generate in each triangle toroid moment directed "up" or "down" by turns and thus build 
up anti-toroic 
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Figure 7b. Triangular anti-regulation of meansquare radii of magnetic moments of each 
domain in (1, 1, 1) plane. 

4 Conclusion 

Now let us remember that there are two electric fields those differ by properties [9]. That, 
which appears in transverse part, is radiation field, unlike the longitudinal part that can 
be connected with the evolution of scalar one in Lorentz gauge 



divA = — ip. 

c 



(4.29) 



In connection with this we should notice the terminological ambiguity. In the framework of 
electrodynamics of continuous media, where may be realized the situation when divj = 0, 
the moments Qi m are the functions, independent to Coulomb moments Qi m , because there 
is no free charges in the system considered. So Qi m cannot be restored as a result of the 
measurement of Coulomb moments by means of permanent electric field. We will illustrate 
it doing inverse transformation, i.e. adding in the current part of multipole expansion the 
transverse and longitudinal contributions, for simplicity Qi m : 

r 1 [Q lm A- 1 Y lm (rt>) - 4 m YUV)A"-MivA], =ln=0 ^ 

=>> QA _1 (-D - VdivA) = QA" 1 (curl curl A - VdivA) = 

c 

= -QA^AA = -QA = QA-^-QA =>■ -QE (4.30) 

Here we take into account that = curl B = curl curl A. Naturally, these expressions 
are nonzero and can be justified only if the system is described but by a charge density 
and Q = / j dr (see e.g. [9]). Note that in quantum mechanics of atoms and moluclues 
there are always isolated charges and the moments Q\ m and Q\ m are connected by the 
evolution equation (3.19) i.e. Qi m = i[H c ,Qi m ]. 
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The Coulomb gauge is generally not applied if a system in consideration contains some 
free charges. Therefore we give the multipole expansion of p(r, t) for the completeness of 
consideration [11, 3] 

= , £ 2^f^TTw /^T o! " i(t)A °^ (r) - (4 ' 31) 

l,m,n=0 V / ' 

<M0 = / ^ +2 "^(r)p(r, t) dr. (4.32) 

Naturally, the Coulomb dipole moments, inherent to (4.24), interact with electric field as 
usual H = — QE. Had we used the Dirac's analysis of constrained Hamiltonian system 
developed within the scope of our problem in [24, 25] we would have been in need of the 
fixation of Coulomb gauge [26] (see also [9]). 
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5 Appendix 

According to definition the current density, corresponding to one dipole moment in the 
proper system of coordinates (see Fig. 1) is equal (see e.g. [13]) 

J Q (q Q , R) = Q a 5{q a - R). (Al) 

where value of dipole is defined as Q Q = — e (q a — R). We may introduce the total 
point-like current J2 a considered to be localized in the fixed point R in the coordinate 
system r and attach it to the fixed point R (Fig. 1): 

J(r) :=£q a *(r-R), (A2) 

a 

where the proper dipoles Q a = — e (q a — R) as before, but Q = q Q . Certainly we 
may obtain the last definition from the exact definition of the total current density of all 
electrons in the system of coordinates r using dipole approximation: 

J(r) = £ cf5(r - q°) « £ q a 5(r - R), (A3) 

a a 

as far as 

5(r-q a ) = 5(r-R + R-q Q )| R _ qa ^5(r - R) + (R - q Q ) W(r - R) + • • • 

But in the system r, the polarization vector field of the molecular system can be repre- 
sented as [14] 

P(r, q Q ) ^£p(q^(r-q Q ) = £ p(q°)5(r-R+R-q a ) | R ^ qa « £ p(q°)5(r-R)+- • • 

a a a 

(A4) 
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where J2 a P a is the sum of dipoles of all electron. From here the polarization current is 
found as 

J(r) :=f =>EP(q°)«(r-R). (A5) 
Comparing (A. 5) with (A. 2) we conclude that 

P(r, q a ) = -e {q a - R) 5(r - R) + • • • (A6) 
Now to define the quadrupole moment we use the following expression 

Qij = \j \PiTj + r,P, - ^(r • P)S ij \ dr. (A7) 
Putting the value P from (A. 6) we obtain 

Qij = -\ [^ r i + r *q? - 2r * r i + jj r (q - r )<%] • ( A8 ) 
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